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Isopiestic vapor pressure measurements have been made on the ternary system water- 
sodium chloride-barium chloride at 25 °C. From these data the activity coefficient of each 
salt in the presence of the other has been calculated. The rule that the logarithm of the 
activity coefficient of one salt (which is equivalent to its chemical potential) varies linearly 
with the concentration of the other salt has been investigated. It is shown that this rule 
is at least approximately valid with the condition usually assumed, namely, that the total 
ionic strength is kept constant. It is shown, moreover, to be valid for another condition, 
namely, that the total concentration of ions is kept constant. This latter condition is more 
suited to isopiestic studies. Finally, the calculation of mutual solubilities from isopiestic 
data is described. 



1. Introduction 

Systems of the type H 2 0-HC1-MC1, where M 
designates a univalent cation, have been studied ex- 
tensively [1, 2], 1 using- concentration cells with hydro- 
gen and silver-silver chloride electrodes. Such cells 
are admirably suited to the determination of the 
chemical potential and activity coefficient of hydro- 
chloric acid in the presence of added chloride and to 
a study of the variation of the activity coefficient of 
hydrochloric acid under selected conditions. Thus, 
one may ask how the activity coefficient of hydro- 
chloric acid varies if the concentration of acid is kept 
constant and sodium chloride is added in increasing 
amounts. Or one may ask how the activity coeffici- 
ent of the acid would vary if both the acid and the 
salt concentration were allowed to change subject to 
a restraining condition; such a condition might be 
that the total ionic strength remained constant. 

Questions of this nature have been answered in 
considerable detail by Harned and his colleagues [1 .]. 
Moreover, studies have been extended to systems in 
which the added chloride is that of a polyvalent 
cation, for example, barium or strontium chloride 
[3, 4] or cadmium chloride [5]. Systems have also 
been studied in which there is no common ion, for 
example, the system water-hydrochloric acid-calcium 
perchlorate [6]. 

When, however, systems containing two salts with 
water as solvent are encountered, this emf method 
has only very limited application. In a few cases, 
suitable electrodes may be available; for example, the 
activity coefficient of zinc chloride in the presence of 
calcium chloride has been measured [7] using a cell 
with zinc amalgam and silver-silver chloride elec- 
trodes. Such well-behaved electrodes are not 
numerous and there are many interesting three- 
component systems where this emf method would be 
difficult, if not impossible, to apply. Thus it is 
difficult to see how suitable electrodes could be 
devised to study the very simple system water- 
sodium chloride-potassium chloride and even more 

1 Figures in brackets indicate the literature references at the end of this paper. 



difficult to suggest an emf method for the system 
water-sodium nitrate-potassium nitrate. 

The isopiestic vapor pressure method has advan- 
tages in this respect and serves as a useful comple- 
ment to the emf method. It has already been used 
to study the systems water-sodium chloride-potas- 
sium chloride [8] and water-potassium chloride- 
potassium bromide [9]. Results are now reported 
for the system water-sodium chloride-barium chlor- 
ide; such results are of interest in themselves because 
they characterize a system of two simple inorganic 
solutes; in addition, this study reveals some difficul- 
ties which have to be overcome in extending the 
isopiestic method to solutions containing salts of 
different valence types. 

2. Methods 

2.1. Definitions 

If an aqueous solution of sodium chloride is in 
isopiestic (vapor-phase) equilibrium with another 
aqueous solution containing both sodium chloride 
and barium chloride, then the two solutions have the 
same aqueous vapor pressure. 

The following definitions are made. 
M B is the molality (moles per kilogram of solvent) of 
a solution containing only sodium chloride (some- 
times termed the "reference" solution). 
M c is the molality of a solution containing barium 
chloride only, in isopiestic equilibrium with the 
sodium chloride solution of molality M B . 
m Bj ni c are the molalities of sodium chloride and 
barium chloride, respectively, in a solution containing 
both these salts and in isopiestic equilibrium with the 
sodium chloride solution of molality M B and the 
barium chloride solution of molality M c . 
R=M B /(rn B +1.5m c ) is the isopiestic ratio. 
#=1.5m c /(m£ + 1.5m c ) is the "ionic fraction" of 
barium chloride in the mixture of sodium chloride 
and barium chloride, m is any linear combination 
of the form 

m=k B m B J rk c m c . 
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For this study it is convenient to put 

k B =l, & c =1.5, m = m B J rl.5m c . 

T B is the (molal) activity coefficient of sodium 
chloride in the reference solution at molality M B in 
equilibrium with the solution containing both 
sodium chloride and barium chloride. 
T c is the activity coefficient of barium chloride at 
molality M c in equilibrium with the reference solu- 
tion of sodium chloride at molality M B . 
7b, Jc are the activity coefficients of sodium chloride 
and barium chloride, respectively, in the mixed solu- 
tion at molalities m B and m c . 

y B is the activity coefficient of sodium chloride in a 
solution containing sodium chloride only at the same 
total ionic strength, I=m B + 3m c , as the mixed solu- 
tion. 7° is the activity coefficient of barium 
chloride in a solution of barium chloride only at the 
same total ionic strength, I=m B +?>m c , as the mixed 
solution. 

<p is the osmotic coefficient of the reference solu- 
tion, defined by the equation 

2M B <p=— 55.51 In a w , 

a w being the water activity of this solution. 

The distinction between T B and y° B and between T c 
and 7° is important. This somewhat complicated, 
but unavoidable, terminology is illustrated schemat- 
ically in figure 1. The primary result of an isopiestic 
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Figure 1. — Schematic representation of the system water- 
sodium chloride-barium chloride. 

The solutions represented by boxes 3, 4, and 5 are in isopiestic equilibrium, 
having the same water activity. The solutions reDresented by boxes 2, 4, and 
6 are not in isopiestic equilibrium but have the same total ionic strength. Mb, 
Mc, yB, Tc are known, as are tub and ra c ; jb and y e , the activity coefficients of 
sodium chloride and barium chloride, respectively, in the solution represented 
by box 4 are to be deduced. 



experiment is the molality M B of a sodium chloride 
solution (box 3) whose water activity a w and activity 
coefficient T B are known. This solution is in vapor- 
phase equilibrium with a barium chloride solution 
(box 5) at molality M c , the activity coefficient of 
barium chloride in this solution, T c , also being known. 
These two solutions are also in vapor-phase equi- 
librium with a mixed solution of sodium chloride and 
barium chloride (box 4) of molalities m B and m c , 
respectively. What are not known are the activity 
coefficients, y B and y c , of sodium chloride and 
barium chloride, respectively, in this mixed solution, 
but it will be shown in the next section that the 
information in these three boxes is sufficient to 
calculate y B and y c . 

One may now enquire how y B is related to y B ; in 
particular, is what is sometimes called Harned's rule 
[1, 2] obeyed? An equation can be written : 



log 7s=log y B —a B I c ~^ B P c 



(1) 



where I c =3 m c and a B , p B , are functions of the total 
ionic strength but not of the individual sodium 
chloride and barium chloride concentrations. 
Harned's Rule is equivalent to the statement that 
/3 B =0. In order to investigate this, a knowledge of 
y° B is needed (box 6) ; that is to say we need to know 
the activity coefficient of a solution of sodium 
chloride only at the same total ionic strength as the 
mixed solution represented by box 4 in figure 1. 
Similarly, for the other component: 



log 7 c =log y c —a c I B —p c F B 



(2) 



where I B =m B , and a knowledge of y° c (box 2) is 
needed to represent the variation of log y c with I B . 

2.2. McKay -Perring Method 

McKay and Perring [10] have, using strict thermo- 
dynamic principles, derived an equation which can 
be put in the form: 



pM B <p 

In y B =\n Y B +\n R+ f(m, M B , x) d (M B <p) 



where 



J v m 1 \d In x/a w m 



m M B 



(3) 

(4) 



The isopiestic ratio, E, is the immediate result of the 
experiment; it is in the evaluation of the integral in 
eq (3) that difficulty is encountered. Let the refer- 
ence solution be equilibrated, not with one mixed 
solution as illustrated in figure 1, but with a series of 
mixed solutions of sodium chloride and barium 
chloride in different proportions, i.e., with different 
values of x, all the solutions having the same value of 
a w . The quantity in the integral, 

1 / dm \ JL(*R\ 

m 2 \d In x)a~ M B \ dx )a w 
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is determined by the tangent to the curve of R versus 
./: for this particular water activity; this tangent will, 
in general, be a function of x and the tangent is to be 
taken at that value of x corresponding to which it is 
desired to evaluate the activity coefficients y B and 
7 C . The experiment must then be repeated to that 
(dR/dx), at this particular value of x, can be deter- 
mined as a function of {M^p) and the integration 
over the interval to (M B <p) performed. There is 
no theoretical way of predicting the variation of 
(d/i/dr) with (M B <p) but a considerable simplification 
ensues if, for the system in question, the isopiestic 
ratio can be expressed as 



R=a—bx—cx 2 



(0) 



where a, 6, and c are functions of M B but not of x\ 
for, in this case, the integral reduces to [8] 



i 



(dM B )d(M B <p) 



and the integral vanishes if c=0, 

R=a—bx. 
In this case, eq (3) takes the very simple form 

\o^y B =lo^T B +logR 
and for the other component 

log7c=logr c +logir 

1.5M C 



where 



(7) 
(8) 
(9) 



R' = 



m#+1.5m c 



2.3. Experimental Procedure 

The apparatus was that used in an earlier investi- 
gation of the system water-sodium chloride-urea 
[11]. Eleven gold-plated silver dishes were used; 
three of them contained the reference sodium chloride 
solution and the remaining eight dishes were divided 
into four groups with four different mixtures each in 
duplicate (or three different mixtures and one pair 
containing a solution of barium chloride only). The 
sodium chloride was part of the sample used pre- 
viously [11]; barium chloride was recrystallized 
twice from water. Analysis showed the presence of 
25 ppm of calcium and 350 ppm of strontium. 2 In 
each run, stock solutions of sodium chloride and 
barium chloride were prepared and mixed in different 
proportions. The stock solution of barium chloride 
was analyzed by the method of Newton and Tip- 
petts [12]. 

The results of one experiment are quoted in detail 
to illustrate the agreement between the molalities 
of solutions in duplicate (or triplicate) dishes : 



2 Analyses by R. W. Burke of the Analysis and Purification Section. 



m B 


2. 6189 
2. 6183 
2. 6182 





1. 7265 
1. 7201 


1. 0075 
1. 0071 


0. 2!) 11 
. 291 1 






m c 


0. 5781 

. 5780 


1. 0450 
1. 0445 


1. 5105 
1. 5107 


1. 6998 
1. 6998 



3. Results 

The experimental results are given in table 1. 
The second column gives the molality (M B ) of the 
reference sodium chloride solution, the third and 

Table 1. Isopiestic data for the system: water — sodium 
chloride (B) — barium chloride (C) at 25 °C 



Set 


Mb 


711 R 


m c 


X 


.R(obs.) 


fl(calc) 


A, 


A 2 


1 


0. 5044 


0. 3591 


0. 1045 


0. 3038 


0. 9780 


0. 9787 


+0. 07 


+0. 06 






.2213 


. 2037 


. 5800 


. 9576 


. 957 1 


-.02 


+.04 






.0794 


. 3060 


. 8520 


. 9308 


. 9304 


-.04 


-.04 








. 3038 


1.0 


. 9244 


.9250 


+.00 







. 5003 


.4033 


. 0736 


0.2149 


.9855 


. 9855 





-.01 






. 2797 


. 1628 


. 4661 


. 9666 


.9662 


-0. 04 


+.01 






.1453 


.2597 


. 7283 


. 9405 


. 9460 


-.05 


-.01 








.3645 


1.0 


.9248 


. 9250 


+.02 


-.04 


2 


. 9728 


0. 73*2 


0. 1644 


0. 2504 


0. 9878 


0. 9880 


+0. 02 









.4910 


. 337H 


. 5079 


.9750 


.9747 


-.03 


-0. 02 






.2610 


. 4999 


.7418 


. 9623 


. 9626 


+.03 


+.04 








. 6832 


1.0 


. 9493 


. 9493 





-.03 


3 


1. 2328 


0. 9479 


0. 1963 


0. 2370 


0. 9923 


0. 992 1 


+0.01 









.6389 


. 4092 


. 4900 


.9841 


. 9S30 


-.05 


-0. 02 






.3356 


. 6193 


.7340 


.9749 


. 9752 


+.03 


+. 05 








.8507 


1.0 


.9661 


.9660 


-.01 


-. 02 


4 


1. 5683 


1. 1653 


0. 2730 


0. 2600 


0. 9959 


0. 9959 





-.01 






0. 7442 


. 5586 


. 5297 


. 9913 


.9913 





+.01 






. 2721 


.8790 


. 8290 


. 9800 


. 9801 


+0.01 


+.01 








1. 0634 


1.0 


. 9832 


. 9832 





-.02 


5 


1. 0226 


1.4141 


0. 1421 


0.1311 


0. 9972 


0. 9978 


+0.00 


+0.08 






0. 9784 


. 4356 


. 4005 


.9944 


. 9930 


-.05 


-.04 






. 5259 


.7414 


. 6790 


. 9900 


.9899 


-.07 


-.07 








1.0983 


1.0 


.9849 


.9853 


+.04 


+.03 


6 


1. 8749 


1. 3930 


0. 3219 


0. 2574 


0. 9995 


0.9991 


-0.04 


-0.01 






1.0035 


.5821 


. 4653 


. 9990 


. 0992 


+.02 


+ . 02 






0. 4308 


. 9647 


.7706 


.9985 


.9994 


+.09 


+.07 








1.2498 


1.0 


1. 0001 


.9995 


-. 00 


-.00 


7 


1. 9873 


1. 2466 


0. 4905 


0.3712 


1. 0025 


1.0027 


+0. 02 


+0.04 






0. 7037 


.8501 


. 6444 


L.0042 


l.ooix 


+. 00 


+ . 00 






.3073 


1. 1124 


. 8445 


1. 0058 


1.0003 


+.05 


+.04 








1. 3149 


1.0 


1.0070 


1.0071 


-.02 


— .04 




2. 0021 


1. 5471 


0. 3000 


0. 2257 


1.0021 


1.0017 


-.04 


—.03 






1. 0496 


. 6303 


. 4739 


1.0035 


1. 0035 





+.01 






0. 5733 


.9452 


. 7121 


1. 0055 


1. 0053 


-0. 02 


— .02 








1. 3245 


1.0 


1. 0077 


1. 0074 


-.03 


— .05 


8 


2. 1724 


1. 8885 


0. 1878 


0. 1298 


1.0010 


1.0011 


+.01 


+.02 






1. 3509 


.5407 


. 3752 


1. 0048 


1.0044 


-.04 


— . 10 






0. 7850 


.9136 


. 6358 


1.0079 


1. 0078 


-.01 


— . 02 






.3184 


1. 2206 


.8519 


1.0107 


1. 0100 


-.01 


+. 06 




2. 1839 


1.0177 


0. 7670 


. 5306 


1. 0072 


1.0004 


-.08 


— . 11 






0. 7602 


. 9375 


. 6491 


1. 0081 


1. 0080 


-.01 


—.02 






. 4557 


1.1377 


. 7892 


1. 0100 


1.0098 


-.02 


+.02 






.3008 


1.2396 


.8607 


1.0110 


1. 0108 


-.02 


+.05 




2. 1980 


1.9192 


0. 1845 


.1260 


1.0010 


1.0011 


+.01 


+.02 






1. 3574 


. 5542 


.3798 


1. 0042 


1. 0044 


+.02 


—.02 






0. 8234 


.9060 


.6227 


1. 0071 


1. 0070 


+0.5 


+.04 






. 3348 


1. 2274 


. 8461 


1. 0102 


1. 0100 


+.04 


+.10 


9 


2. 3615 


1. 7350 


0. 4108 


0. 2621 


1. 0044 


1.0043 


-0. 01 









1. 1799 


.7744 


. 4961 


1. 0085 


1. 0080 


+.01 









0. 5922 


1. 1591 


.7459 


1. 0131 


1.0131 














1. 5466 


1.0 


1. 0179 


1. 0178 


-0.01 





10 


2. 6185 


1. 7263 


0. 5781 


0. 3343 


1. 0097 


1.0098 


+0.01 









1.0073 


1.0448 


. 6087 


1.0171 


1. 0170 


-.01 









0.2911 


1. 5100 


. 8862 


1. 0241 


1. 0241 





+0.01 








1. 6998 


1.0 


1. 0270 


1. 0271 


+0.01 


— . 01 


11 


2. 7147 


2. 0250 


0. 4451 


0. 2480 


1. 0082 


1.0087 


+.05 









1. 3969 


.8514 


.4776 


1. 0152 


1.0147 


-.05 


—0.01 






0. 6356 


1. 3458 


. 7605 


1. 0228 


1. 0221 


-.07 


—.03 








1. 7610 


1.0 


1. 0277 


1. 0283 


+.00 


+.03 



Ai = ii?(calcbyeq(6))-7^(obs)lX102. 
A 2 =[i?(calc by eq(5))-i?(obs)]X102 
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fourth columns give values of the molalities of sodium 
chloride (m B ) and barium chloride (m c ) in the mixed 
solutions which were found to be in isopiestic 
equilibrium with the reference sodium chloride 
solution. The next column gives values of the 
"ionic fraction' ' (x) and the sixth values of the 
isopiestic ratio. Figure 2 is a plot of isopiestic ratio 
versus ionic fraction for 3 of the 11 sets of results 
obtained; from this plot it would appear that R is a 
linear function of x. The data in each set were 
therefore fitted to eq (7) and the a and b parameters 
evaluated by the method of least squares. These 
values of a and b are given in the second and third 
columns of table 2 and the values of R calculated by 
eq (7) are given in the seventh column of table 1. 
The next column gives the difference (multiplied by 
100) between the observed values of R and those 
calculated by eq (7). The average of these differ- 
ences for each set is given in the sixth column of 
table 2. It will be seen from figure 2 and from table 
2 that the a parameter does not have exactly the 
theoretical value of unity, suggesting that the plot 
of R versus x does have slight curvature, more 
pronounced at low values of x. The data were 
therefore fitted to eq (6) with a=l and the b and c 
parameters, calculated by the method of least 
squares, are given in table 2. Values of R have been 
calculated with the aid of these parameters; they are 
not recorded in table 1 but the last column of this 
table gives the difference (multiplied by 100) be- 
tween the observed values of R and those calculated 
by eq (6) ; the average difference for each set is given 
in the last column of table 2. It will be seen that 
any superiority eq (6) may have over eq (7) in repre- 
senting the experimental data is within the experi- 
mental error; the simpler equation, (7), will therefore 
be used in subsequent calculations. 



I.UD 


! 1 

^^--O SET II 




^^^-cr^""^ —-0 _ __ __ O SET 7 


1.0 




0.95 


X. 




>0 SET 1 


0.9 


1 1 



Table 2. Values of the parameters of eq (7) and eq (6) 





Eq(7) 


Eq(6 


, a=l 






Set 










Avg 

Ai 


Avg 












A 2 




a 


b 


6 


c 






1 


1. 0021 


0. 0771 


0. 0681 


0. 0075 


0.04 


0.03 


2 


1. 0009 


.0516 


.0481 


. 0029 


.02 


.02 


3 


1.0006 


.0345 


.0318 


. 0023 


.03 


.03 


4 


1.0004 


.0172 


.0156 


.0014 





.01 


5 


0. 9997 


.0144 


.0152 


-. 0004 


0.06 


.05 


6 


.9990 


.0006 


-. 0030 


. 0025 


.05 


.04 


7 


1. 0000 


-.0074 


-. 0080 


.0008 


.03 


.04 


8 


0. 9994 


-.0132 


-. 0085 


-. 0056 


.04 


.05 


9 


.9995 


-.0183 


-. 0164 


-. 0015 


.01 





10 


1.0012 


-. 0259 


-. 0300 


.0031 


.01 


0.01 


11 


1.0023 


-.0260 


-. 0349 


.0069 


.06 


.02 



4. Calculation of Activity Coefficients 

This study of a mixture of a 1:1 electrolyte and a 
2:1 electrolyte reveals a difficulty not encountered 
with mixtures of two electrolytes of the same 
valence type. This can be illustrated as follows. 
It was found [8] in a study of the system water- 
sodium chloride-potassium chloride that the follow- 
ing solutions were in isopiestic equilibrium: 





^NaCI 


2. 


3248 


1. 


8865 


1. 


4490 


1. 


0324 


0. 


5630 




^KCl 







0. 


4894 


0. 


9701 


1. 


4238 


1. 


9269 


%aC 


i+mKCi 


2. 


3248 


2. 


3759 


2. 


4191 


2. 


4562 


2. 


4899 



Thus only a small interpolation was needed to study 
the system at a total ionic strength of 2.4. It is 
very different with the system water-sodium chloride- 
barium chloride as the following example, taken from 
set 9 of table 1, shows: 



m B 

m c 

ras + 3 m c 

%+ 1.5 m c 



2. 3615 


1. 7350 


1. 1799 


0. 5922 





0. 4108 


0. 7744 


1. 1591 


2. 3615 


2. 9574 


3. 5031 


4. 0695 


2. 3615 


2. 3512 


2. 3415 


2. 3309 





1. 5466 
4. 6398 

2. 3199 



Thus what remains approximately constant is not 
the total ionic strength (m B +3m c ) but the "total 
ionic concentration " as measured by (m fi +1.5 m c ); 
for, as the above table shows, the total ionic con- 
centrations of the five isopiestic solutions are very 
nearly the same — the average value is 2.3410 and 
the largest deviation from this value is only 0.8 per- 
cent. We therefore considered whether Harned's 
Rule can hold in this form: the logarithm of the 
activity coefficient of sodium chloride is a linear 
function of the molality of barium chloride in mix- 
tures of constant total ionic concentration, with a 
corresponding proposition for the activity coefficient 
of barium chloride. We therefore rewrite eq (1) and 
eq (2) as: 



log y B =log y B —a B m c 



(10) 
(11) 
Figure 2. Isopiestic ratio (R) versus ionic fraction (x). ' 7°b and y° c do not have the meanings assigned to 



log 



c =lo£ 



7 C — a c m B 
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them in section 2.1; in the present context, y B is now 
the activity coefficient of sodium chloride in a solu- 
tion of sodium chloride only at a molality equal to 
the ionic concentration (m 7i +1.5 m c ) of the mixed 
solution; 7° that of barium chloride in a solution of 
barium chloride only at a molality equal to (m s +1.5 
m c )/1.5. It would be possible to devise new symbols 
but there is no need to complicate the symbolism 
still further. 

We now use eq (8) and eq (9) to calculate the 
activity coefficients of each salt in the mixture at 
constant total ionic concentration and then calculate 
a B and a c by eq (10) and eq (11). Values of y B have 
already been tabulated [13]; data are available for 
barium chloride but the osmotic and activity co- 
efficients have been computed again 3 from earlier 
isopiestic data (14,15] together with those results in 
table 1 which pertain to solutions of barium chloride 
without added sodium chloride. The data are given 
in table 3. 



Table 3. Osmotic and activity coefficients of barium chloride 
in aqueous solution at 25 °C 



m 


<fc 


Te 


m 


<Pc 


7c 


0.10 


0. 8422 


0. 5078 


1.10 


0. 9520 


0. 4076 


.20 


. 8382 


. 4523 


1.20 


.9681 


.4128 


.30 


.8440 


. 4265 


1.30 


.9845 


.4188 


.40 


.8534 


.4121 


1.40 


1.0010 


.4255 


.50 


.8648 


.4039 


1.50 


1. 0175 


. 4329 


.60 


.8775 


.3996 


1.60 


1.0340 


.4408 


.70 


. 8912 


.3979 


1.70 


1. 0503 


. 4492 


.80 


.9056 


. 3982 


1.80 


1. 0665 


.4580 


.90 


. 9206 


.4001 


1.90 


1. 0823 


. 4672 


1.00 


.9361 


.4033 


2.00 


1. 0978 


.4767 



Table 4 gives the values of a B and a c . Sub- 
stituted in eq (10) and eq (11), they will represent 
the data with a maximum difference of 0.0005 in 
log y B or log y c , the average difference being much 
less than this. Within experimental error, there- 
fore, it can be concluded that the proposed modifica- 
tion of Harned's Rule is valid for this system. Figure 
3 is a plot of a B and a c versus M B . This emphasizes 
the fact that a B and a c have large values in dilute 
solutions and small values in concentrated solutions. 
Thus there are large excess free energy terms in mix- 
ing dilute solutions of sodium chloride and barium 
chloride at constant total ionic concentration but 
these terms are much smaller (and of opposite sign) 
if the solutions are concentrated. This is a some- 
what surprising result but confirms what has been 
observed before, particularly for the system water' 
hydrochloric acid-cesium chloride. 

Table 4 also contains values of (a 5 +1.5 a c ) and 
these are plotted versus M B in figure 3. Now it has 
been shown [16, see also appendix] as a direct 
result of the cross-differentiation relation, 

/MogrA =3 / Qlog7A (12) 

V &m c / mn \ dm B J mr 



3 We thank Dr. Wallace Davis, Jr., of the Oak Ridge National Laboratory for 
these computations. 



that the a B and a c parameters of eq (10) and eq (11) 
must obey the relation 

(a B -\- 1.5 a f )=constant independent of (7^ + 1.5 m c ) 

(13) 
provided that 

(1) a B and a c are functions of (rn B -\- 1.5 m c ) 
only, and 

(2) there are no terms in m B and m c other than 
a B m c and a c m B , i.e., terms such as @ B m 2 c are excluded. 
It will be seen that (a B -\-l.5 a c ) is very nearly 
constant, the largest deviation is at the lowest 
concentration where (a B + 1.5 a c ) is the summation 
of two large terms of opposite sign. No significance 
is attached to the very slight variation of (a 7i + 1.5 a c ) 
which is well within the experimental error and it 
is concluded that eq (10) and eq (11) suffice without 
the need for additional terms in higher powers of 
m B and m c . 

Table 4. Values of the parameters an and a r of eg (10) and 
eq (11) 



Set 


aB 


OLc 


(ccB-r-1.5a e ) 


1 


0. 0859 


-0. 0602 


-0. 0044 


2 


. 0320 


- .0225 


- .0018 


3 


.0175 


- .0128 


- .0017 


4 


. 0069 


- .0047 


- .0002 


5 


. 0060 


- .0046 


- .0000 


6 


. 0005 





. 0005 


7 


-. 0025 


. 0019 


.0004 


8 


-.0040 


. 0033 


.0010 


9 


-. 0055 


.0044 


.0011 


10 


-. 0081 


0058 


.0000 


11 


-. 0085 


.0055 


- .0002 




Figure 3. a B , a Cf and (a B -\-1.5a e ) versus M B for systems at 
constant total ionic concentration. 
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Figure 4. Span of total ionic strength in each set of 
measurements. 



We could stop at this point because we have 
obtained a complete description of the system over 
the concentration range 0.5 <^(m 5 +1.5 m c )<^ 2.7; 
eq (10) and eq (11) with the a B , a c parameters of 
table 4 give this information. It is of interest, 
however, to pursue the calculation and consider 
variations at constant total ionic strength instead 
of constant total ionic concentration. The difficulty 
of applying the isopiestic method to systems of 
different valence types at constant total ionic 
strength has already been referred to; figure 4 
illustrates this in another way by showing the span 
of values of total ionic strength covered in each of 
the eleven sets of measurements. It shows that, 
whereas the set of five results quoted above for the 
system water-sodium chloride-potassium chloride 
contribute five points to the data at a total ionic 
strength, 2.4, any one set in the water-sodium 
chloride-barium chloride system contributes ef- 
fectively only one point at any selected total ionic 
strength. Thus, results from different sets have to 
be harnessed in order to obtain data at different 
values of x at any given total ionic strength. More- 
over, the number of contributions depends on the 
total ionic strength. Thus, at 7=3, sets 4 to 11 
make eight contributions; at 1=5, sets 10 and 11 
only contribute. 

We proceed as follows. At a total ionic strength 
of 7—3, for example, we could have five solutions of 
the following compositions : 



No. 


m B = I B 


Ic 


1.5 m c 


m 


1 


0. 5 


2.5 


1. 25 


1.75 


2 


1.0 


2.0 


1.00 


2.00 


3 


1. 5 


1. 5 


0. 75 


2. 25 


4 


2.0 


1.0 


.50 


2.50 


5 


2. 5 


0. 5 


.25 


2.75 



A value of log y B for each solution can be calculated 
by eq (10) using the value of log y° B appropriate to 
each value of m and a value of a B for each value of ra 
is read from a large scale plot of figure 3. Then 
a new value of a B is calculated by eq (1), ignoring 
the fi B P c terms. The calculation is as follows: 



No. 


-log 7b 


a B eq (10) 


-logTs 


a B eq (1) 


— log 7 s 

(calc.) 


1 


0. 1799 


0. 0025 


0. 1820 


0. 0142 


0. 1818 


2 


. 1755 


-. 0018 


. 1743 


. 0139 


. 1747 


3 


. 1702 


-. 0050 


. 1677 


. 0141 


. 1677 


4 


. 1627 


-. 0070 


. 1604 


.0139 


. 1606 


5 


. 1550 


-. 0085 


. 1531 


.0132 


. 1536 



The best value of a B for use in eq (1) is 0.0141 and, 
substituting for —log 7^ = 0.1465 at 7=3, we get 

log7 B =— 0.1465— 0.0141 I c 

Values of —log y B calculated by this equation are 
given in the last column of the above table; it is 
evident that, at least at 7=3, eq (1) is adequate 
with only a term linear in I c , so that /3 B = 0. Table 5 

Table 5. Values of the parameters of eq (1) and eq {2) 





0b= 


/3e=0 




I 


OLB 


Ctc 


(2as+<*c) 


0. 75 


0. 0082 


-0. 0077 


0. 0087 


1.5 


. 0095 


-.0085 


. 0105 


2.0 


.0111 


-.0118 


.0104 


2.5 


. 0127 


-.0155 


.0099 


3.0 


. 0141 


-.0173 


. 0109 


3.5 


.0159 


-.0191 


. 0127 


4.0 


.0167 


-.0206 


.0128 


4.5 


.0178 


-.0218 


.0138 


5.0 


.0188 


-.0225 


.0151 



gives values of a B and a c calculated in this way at 
different values of the total ionic strength. The last 
column gives values of (2a B -\-a c ), a quantity which 
is important because it can be shown [16, see also 
appendix] that in a system represented by eq (1) 
and eq (2) with (3 B =p c =0, the quantity (2a B J ra c ) 
should be independent of the total ionic strength. 
Not all the values of a B and a c in table 5 are of equal 
reliability, partly because at some values of / we 
can use more sets of data from tables 1 and 4 than 
at others, and partly because there is some difficulty 
in interpolating in figure 5 at low values of /; the data 
at 1=2.5, 3.0, and 3.5 are probably more reliable, 
those at higher and lower total ionic strengths less so. 
Figure 5 compares these values of a with those 
given by Harned and Gary [1, 4] for the system water- 
hydrochloric acid-barium chloride. We note that the 
a value for barium chloride is negative in both sys- 
tems, so that the replacement of barium chloride by 
either hydrochloric acid or sodium chloride in a 
system at constant total ionic strength raises the 
activity coefficient of barium chloride; the effect is, 
however, much larger with hydrochloric acid. The 
a values for hydrochloric acid and sodium chloride 
are positive so that the activity coefficient of each is 
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+ 0.1 



+ 0.05 — 




-0.05 — 



Figure 5. Parameters of eq (1) and eq {2). 

System water; sodium chloride (B)-barium chloride (C) curve 1 as; curve 
2a r "; curve 3 (2aj»+a e ), the curve lias been raised by 0.01 for the sake of clarity. 

System water-hydrochloric acid (B)-harium chloride (C) curve 4 a b; curve 
5 a c ; curve 6 (2 aB+a c ). 



decreased by barium chloride; again the effect is 
very much larger with hydrochloric acid. For the 
sodium chloride-barium chloride system, (2a B J r a c ) 
increases slightly with increasing total ionic strength. 
The effect may be spurious, due to the difficulty we 
have in evaluating data in systems of constant total 
ionic strength. If the effect is real, it suggests that 
the /3 terms in eq (1) and eq (2) are not both zero; 
B =— 0.0003, C =O would satisfy the condition 
demanded by the cross-differentiation relation, eq 
(12). 

5. Solubility Relations 

Let m s be the solubility of sodium chloride in 
water and 7, its activity coefficient at this molality. 



Let m B be the solubility of sodium chloride in a 
solution of barium chloride of molality m c and y B 
its activity coefficient. Then 



fn 



or 
or 



m 2 s y 2 =m B (m B +2m c )y 2 1 

m s y s =m B (l+2m c /m B y /2 y B 

log m B /m,= — log 7b/t»— 1/2 log (l+2m c /m B ). 

Furthermore, 

log Ts^log 7°b— 3a B m c . 

At 25 °C, m 8 = 6.147 mole kg -1 [17] and log y s = 
0.0020. Moreover, in the vicinity of m B = 6 mole 
kg" 1 , log y° B is almost linear in the total ionic 
strength, 

log 7i=log 7.S+0.0545 (m B +3m c — m 8 ). 

Hence 

log m B /m, s = — 0.0545 {m B +3m c — m s )+3a B m c 

-1/2 log (l+2m e /m B ). 

Now a B is a function of (m B +3m c ) and m B is not 
known in advance. However, an inspection of 
table 5 and figure 5 shows that in concentrated 
solutions the value of a B is not very sensitive to 
the value of the total ionic strength and for our 
present purpose we can put a^ = 0.02 as a good ap- 
proximation. The solution of this equation predicts 
a solubility of 5.465 mole kg" 1 in 0.5 w barium 
chloride solution and of 4. SI 3 mole kg l in 1.0 m 
barium chloride solution (such solution would be 
supersaturated with respect to barium chloride). 

The solubility of barium chloride in sodium 
chloride solutions is governed by somewhat different 
considerations because it is the dihydrate which is 
the stable solid phase. Representing the molality 
and activity coefficient of barium chloride in the 
saturated aqueous solution by m, and 7,, respec- 
tively, we have : 

4m 3 s y*al {s) = m c (m B +2m c ) 2 y 3 a 2 w 

where a w ( 8 ) is the water activity of the saturated 
aqueous solution and a w that of the saturated solu- 
tion containing sodium chloride at molality m B . 
Thus, 



log m c /m 8 =— logy c /y s — 2/3 log (l + m B /2m f ) 

—2/3 log (a w /a W ( S) ). 

At 25 °C,m s = 1.790 mole kg^tlS], log 7,= -0.3400. 
Again, log 7° is almost linear in the total ionic 
strength : 

log 7c=log Y a +0.0283 (m B +3m c — 3m s ) 
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so that 

log m c /m s =— 0.0283 (m B + 3m c — 3m s )+a c m B 

-§ log (l+m B /2 m c ) — § log (aja wis) ). 

From table 5 and figure 5 it would appear that a c ~ — 
0.023 for the saturated solution. In this way we 
have calculated solubilities of barium chloride of 
1.401, 1.045, 0.736, 0.484, 0.299, 0.175, and 0.100 
moles/kg in solutions of 1.0, 2.0, 3.0, 4.0, 5.0, 6.0, 
and 7.0 r/i-NaCl, respectively. These data have 
been used to construct the two curves in figure 6. 
There seem to be no direct measurements of mutual 
solubility at 25° but Findlay and Cruickshank [19] 
determined the compositions of eight saturated so- 
lutions at 20 °C; these are represented by closed 
circles in figure 6. Schreinemakers and de Baat 
[20] made eight measurements at 30 °C and these are 
represented by open circles in figure 6. In spite of 
the approximations we made in calculating solubili- 
ties from isopiestic data, the agreement is good and 
confirms the validity of the data given in table 5 to 
describe the activity coefficients of the salts in the 
system water-sodium chloride-barium chloride. 



6.0 



5.0 



4.0 



SOLUBILITY OF NaCl 
IN B0Cl 2 SOLUTIONS 




3.0 



2.0 



SOLUBILITY OF BOCl 2 
IN NQCl SOLUTIONS " 



Figure 6. 



m r 



Solubility relations in the system water-sodium 
chloride-barium chloride. 



, Solubility curves calculated from isopiestic data at 25 °C. 

#, Direct measurements at 20 °C. 
0» Direct measurements at 30 °C. 



7. Appendix 

Although the derivation of eq (13) has already 
been described [1, 16], it is worthwhile considering 
the more general case of a mixture of electrolytes of 
different valence types for which the terms of eq 
(1) and eq (2) are not zero. Consider a solution of 
electrolyte B (giving v B ions) at molality m B , and 
electrolyte C (giving v c ions) at molality m c . Let 
f be a linear combination of the form 

{=k B m B +k c m c 

and let the activity coefficients at constant f be 
given by 

log Ys=log y° B —A B m c —B B m 2 c 

log Yc=log y c —A c m B —B c m B 

where A B and A c are functions of f but not of m B 
and m c individually and B B , B c are parameters inde- 
pendent of f , m B , and m c . The cross-differentiation 
relation : 



v B 



/ d log y B \ =v / dlogTA 
\ dm c J mn Vc \ dm B J n 



can be written: 

/ d log y B \ _ , / d log y c \ 

Moreover, 

log 75=log y° B —A B (£—k B m B )/k c —B B ({—k B m B ) 2 /k 2 c 
and 

log 7c=log y° c — A c ({— k c m c )/k B — 5 c (f — k c rn c ) 2 /k 2 B . 
As log y° B , log 7c, A B) and A c are functions of f only, 



V B i 



*{ 



d log y° B A B t—k B m B dA B 



= v c k. 



d{ k c 

d log y° c 



h 



■{ 



A e 



d£ k B 



-TUi 



dA c 



-2B B 



t-k 



-JCsMb \ 

k\ J 



-2B r m 



s /k B V- 



If this equality is to hold for all values of f and m B , 
the coefficient of m B must vanish, and 



k B k t 



V B K B > 



or 



B -\-v c k B k c 



# 



dAc 



= — 2(v B k B B B +v c k c B c ) 

k B k c {v B A B + v C A c ) = constant 
— 2(v B k B B B + v c k c B c ){ 



If v B = 2, u c =3, k B =l, k c =3, as with sodium chloride 
and barium chloride at constant total ionic strength, 
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(6A B + 9^1 C ) = constant — 2(2B B + 9B C )I or, putting 
A B = \\a Bl A c =a c , B B = 9(3 B , B C = /3 C , 

Log y B =\og y° B —a B I, — fi B I 2 c 
loo- y c = log y° c — a c I B —(5 c P B 

which are eq (1) and eq (2), then 

(2a B +a c ) = constant - 2 (2ft, + C ) 7. 

If 1^=2, u c =3, £# = 1, & c =1.5, as with sodium 
chloride and barium chloride at constant total ionic 
concentration, 

1.5 {2 An + 3^4 c ) = constant -2(2^+4.5£ c ) (m B + 
1.5m c ) 

or, putting ^4^ = ^, A c =a CJ t3 B = p c =Q, 

log75=log y°B—a B m c , 
logy c =logy° c —a c m B 

which are eq (10) and eq (11), then 

(<x B -\-1.5a c ) = constant. 

Alternatively, we could put A B =1.5 a B , 

A c =a c , B B =\.5 2 I3 BJ B c = p c , so that: 

logYB=log7s— a B (1.5m e )— fe(1.5m c ) 2 ? 

log T c =log y° c — a c m B — p c m 2 B 

and obtain the more compact equation : 

(a B +a c ) = consUmt — 2 (P B + p c ) (m B + 1.5m c ) . 
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